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Abstract. The paper is concerned with patchy vector fields, a class of discontinuous, piecewise 
smooth vector fields that were introduced in [A-B] to study feedback stabilization problems. We 
prove the stability of the corresponding solution set w.r.t. a wide class of impulsive perturbations. 
These results yield the robusteness of patchy feedback controls in the presence of measurement errors 
and external disturbances. 
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1 - Introduction and Basic Notations. 

Aim of this paper is to establish the stability of the set of trajectories of a patchy vector field 
w.r.t. various types of perturbations, and the robustness of patchy feedback controls. 

Patchy vector fields were introduced in [A-B] in order to study feedback stabilization problems. 
The underlying motivation is the following: The analysis of stabilization problems by means of 
Lyapunov functions usually leads to stabilizing feedbacks with a wild set of discontinuities. On the 
other hand, as shown in [A-B], by patching together open loop controls one can always construct 
a piecewise smooth stabilizing feedback whose discontinuities have a very simple structure. In 
particular, one can develop the whole theory by studying the corresponding discontinuous O.D.E's 
within the classical framework of Caratheodory solutions. We recall here the main definitions: 

Definition 1.1. By a patch we mean a pair (f2, g) where VL C W 1 is an open domain with smooth 
boundary dil, and g is a smooth vector field defined on a neighborhood of the closure ft, which 
points strictly inward at each boundary point x € dfl. 

Calling n(x) the outer normal at the boundary point x, we thus require 

(g{x), n(x)) < for all x E dQ. (1.1) 



Definition 1.2. We say that g : O i— > M. n is a patchy vector field on the open domain if there 
exists a family of patches g a ); a E A} such that 

- A is a totally ordered set of indices; 

- the open sets fi Q form a locally finite covering of f2, i.e. £1 = U aeJ {Q a and every compact set K C W 1 
intersect only a finite number of domains Q a , a E A; 

- the vector field g can be written in the form 

g(x)=g a {x) if xen a \\jQp. (1.2) 

I3>a 

By setting 

a*(x) = max {a E A ; x E £1 Q }, (1.3) 
we can write (1.2) in the equivalent form 

g(x) — g a , (x) (x) for all ie!l. (1.4) 

We shall occasionally adopt the longer notation (f2, g, (Cl a , g a ) aeA ) to indicate a patchy vector 
field, specifying both the domain and the single patches. If g is a patchy vector field, the differential 
equation 

x = g(x) (1.5) 

has many interesting properties. In particular, in [A-B] it was proved that the set of Caratheodory 
solutions of (1.5) is closed in the topology of uniform convergence, but possibly not connected. 
Moreover, given an initial condition 

x(t )=x , (1.6) 
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the Cauchy problem (1.5)-(1.6) has at least one forward solution, and at most one backward solution. 
For every Caratheodory solution x = x{t) of (1.5), the map t i— > a*(x(t)) is left continuous and 
non-decreasing. 

In this paper we study the stability of the solution set for (1.5) w.r.t. various perturbations. 
Most of our analysis will be concerned with impulsive perturbations, described by 

y = g(y) + w. (1.7) 

Here w = w(t) is any left continuous function with bounded variation. By a solution of the perturbed 
system (1.7) with an initial condition 

y(to) = 2/o, (1.8) 
we mean a measurable function t \— > y(t) such that 

y(t) = y + [ g(y(s)) ds + [w(t) - w(t )] . (1.9) 

J to 

If w(-) is discontinuous, the system (1.7) has impulsive behavior and the solution y(-) will be discon- 
tinuous as well. We choose to work with (1.7) because it provides a simple and general framework 
to study robustness properties. Indeed, consider a system with both inner and outer perturbations, 
of the form 

x = g(x + e 1 (t))+e 2 (t). (1.10) 
The map t \— ► y(t) = x(t) + e\{t) then satisfies the impulsive equation 

y = g(y) + e 2 (i) + ei(i) = g(y) + w, 
where t 

w(t) = ei(t) + / e 2 (s)ds. 

J t 

Therefore, from the stability of solutions of (1.7) w. r. t. small BV perturbations w, one can im- 
mediately deduce a result on the stability of solutions of (1.10), when Tot.Var.jei} and I^Hl 1 
are suitably small. Here, Tot.Var.jei} denotes the total variation of the function e\ over the 
whole interval where it is defined, while Tot.Var.jei ; J} denotes the total variation of e\ over 
a subset J. Any BV function w — w(t) can be redefined up to L 1 -equivalence. For sake of def- 
initcness, throughout the paper we shall always consider left continuous representatives, so that 
w(t) = w(t~) = lim w(s) for every t. The Lebesgue measure of a Borel set J C M will be denoted 

s— >t — 

by meas(J). 

We observe that, since the Cauchy problem for (1.5) does not have forward uniqueness and 
continuous dependence, one clearly cannot expect that a single solution of (1.5) be stable under 
small perturbations. What we prove is a different stability property, involving not a single trajectory 
but the whole solution set: If the perturbation w is small in the BV norm, then every solution of 
(1.7) is close to some solution of (1.5). This is essentially an upper scmicontinuity property of the 
solution set. 



Theorem 1. Let g be a patchy vector field on an open domain !1 Cl". Consider a sequence of 
solutions y u { ) of the perturbed system 

yu = g(y„) + w u ie[0,T], (1.11) 
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with Tot.Var.{w„} — ► as v — > oo. If the y v : [0,T] converge to a function y : [0,T] i— > O, 
uniformly on [0,T], then y(-) is a Caratheodory solution of (1.5) on [0, T]. 

Corollary 1.3. Lef g be a patchy vector field on an open domain f2 C W l . Given any closed subset 
icfi, any compact set K C A, and any T, e > 0, there exists S = 5{A, K,T,e) > such that the 
following holds. If y : [0,T] ^ A is a solution of the perturbed system (1.7), with y(0) G K and 
Tot.Var.{w} < S, then there exists a solution x : [0,T] ^ f2 of the unperturbed equation (1.5) with 

\\ X -vL°o ([0 ,T\) <£ ( L12 ) 

We remark that the type of stability described above is precisely what is needed in many 
applications to feedback control. As an example, consider the problem of stabilizing to the origin 
the control system 

x = f(x,u). (1.13) 

Given a compact set K and e > 0, assume that there exists a piecewise constant feedback u — U(x) 
such that g(x) = f{x, U(x)) is a patchy vector field, and such that every solution of (1.5) starting 
from a point x(0) G K is steered inside the ball B e centered at the origin with radius s, within a time 
T > 0. By Corollary 1, if the perturbation w is sufficiently small (in the BV norm), every solution 
of the perturbed system (1.7) will be steered inside the ball B 2e within time T. In other words, the 
feedback still performs well in the presence of small perturbations. Applications to feedback control 
will be discussed in more detail in Section 3. 



2 - Stability of Patchy Vector Fields. 

We begin by proving a local existence result for solutions of the perturbed system (1.7). 

Proposition 2.1. Let g be a patchy vector field on an open domain C K™. Given any compact 
set K C fl, there exists \ = \k > such that, for each yo € K, to G M, and for every Lipschitz 
continuous function w — w(t), with Lipschitz constant ||w||l°° < X, the Cauchy problem (1.7)-(1.8) 
has at least one local forward solution. 

Proof. 

Fix some compact subset K' C O whose interior contains K . To prove the local existence of a 
forward solution to (1.7), first observe that, because of the inward-pointing condition condition (1.1) 
and the smoothness assumptions on the vector fields g a , one can find for any a G A some constant 
Xa > such that 

sup (g a (x) + v, n a (x)) < 0, (2.1) 

\v\<Xa 

where n a (x) is the outer normal to dil a at the boundary point x. Since K' is a compact set and 
{fl a }a is a locally finite covering of fi, there will be only finitely many elements of {f2 a } a that 
intersect K' . Let 

{ai,.. .,a N } = {aeA : Sl a r\K'?Q}, (2.2) 
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and, by possibly renaming the indices c^, assume that 

ai <•••<£*„. (2.3) 

Choose a constant \ > such that 

A<inf{xa s : i = l,...,N}. (2.4) 
For any fixed yo E K, consider the index 

S(j/ ) = max {a : y E fl a } ■ 
By the definition of x, any solution y = y(-) to the Cauchy problem 

y = 9a(y) + w, y{to) = yo, 

associated to a piecewise Lipschitz map w = w(t) with ||to||i,°° < X, remains inside Qa for all 
t E [to, to + 5], for some S > 0. Hence, it provides also a solution to (1.6) on some interval [to, to + S'], 
< 5' < 5. □ 



Toward a proof of Theorem 1, we first derive an intermediate result. By the basic properties of a 
patchy vector field, for every solution t \— > x(t) of (1.5) the corresponding map t \— > a*(x(t)) in (1.3) 
is nondecreasing. Roughly speaking, a trajectory can move from a patch £l a to another patch Clp 
only if a < j3. This property no longer holds in the presence of an impulsive perturbation. However, 
the next proposition shows that for a solution y of (1.7) the corresponding map 1 1— > a* (y(t)) is still 
nondecreasing, after a possible modification on a small set of times. Alternatively, one can slightly 
modify the impulsive perturbation w, say replacing it by another perturbation , such that the 
map 1 1— > a*(y^(t)) is monotone along the corresponding trajectory t \— > j/^(t). 

Proposition 2.2. Let g be a patchy vector field on an open domain ft C M™, determined by the 
family of patches {(Q a , <? Q ); a £ .4}. For any T > and any compact set K C fi, t/iere exist 
constants C,8 > and an integer N such that the following holds. 

(i) For every w E BV wzi/i Toi. Var.{tu} < 5, and /or eijery solution y : [0, T] il of the Cauchy 
problem (1.7)-(1.8) with y E K, there is a partition of [0, T], = T\ < t 2 < ■ ■ ■ < tn+i = T, 
and indices 

a\ < a 2 < ■ ■ ■ < ol 7v , (2.5) 

such that 

a*(y(t))>a t V t e]n, n+i], i = 0,...,N, (2.6) 

meas ( |J {t G [n, T i+1 ] : a*{y{t)) > a,}) < C • Tot.Var.{w}. (2.7) 

i>0 



(ii) For every BV function w = w(t) with Tot.Var.{w} < S, and for every solution y : [0,T] i— » Ct 
of the Cauchy problem (1.7)-(1.8) with y E K, there is a BV function = tu^(t) and a 
solution y^ : [0, T] ^ f2 of 

y^=g(y^)+w^, (2.8) 
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so that the map t \— > a*(y^(i)) is non- decreasing and left continuous, and there holds 

Tot. Var.iw^} < C ■ Tot. Var.{w} , 



(2.9) 

<T n . Tn+ \7nv fn„\ 

L~([0,T]) 



l^-fllrocrro™ Tot.Var.{w}. 



Proof. 

1. The proof of (i) will be given in three steps. 

Step 1. Since each g a is a smooth vector field and we are assuming a uniform bound on the 
total variation of every perturbation w = w(t), there will be some compact subset K' C fi that 
contains every solution y : [0,T] i— ► Q of (1.7) starting at a point yo € K. We will assume without 
loss of generality that every domain Q, a is bounded since, otherwise, one can replace fl a with its 
intersection Q a n 0' with a bounded domain (1'cfi that contains K , preserving the inward-point 
condition (1.1). For each aei, define the map <p a : i— ► R by setting 

{d(x, dil a ) if ieH a , 

(2.10) 
— d(x, dfl a ) otherwise, 

and let 

ip+ (x) = max{^ a (i), 0} 

denote the positive part of ip a (x). The regularity assumptions on the patch £l a guarantee that 
(p a is smooth if restricted to a sufficiently small neighborhood of the boundary dVL a . Thus, if 
: i = l,..., N} denotes the finite collection of domains that intersect the compact set K' as 
in (2.2)-(2.3), there will be some constant ~p > so that, setting 

nl = {xen ■. d{x, on a ) >p}, (2.11) 

the restriction of any map ip ai to the domain \ fig,, be smooth. In particular, for any i = 1, . . . , N, 
we will have 

Vy ai {x) = -n ai (7r ai (a;)), Vie(]\^, (2.12) 

where n a . represents as usual the outer normal to dtt ai , while 7r Qi (x) denotes the projection of the 
point x onto the set dfl ai . On the other hand, thanks to the inward-pointing condition (1.1), we can 
choose the constant p so that 

sup (g ai {x), n ai (n ai (x))) < -c', (2.13) 

i=l,...,JV 

for some d > 0. Moreover, the smoothness of the fields g a on £1 implies the existence of some c" > 
such that 

sup (g aj {x), n a .(7r Qi (x))) < c". (2.14) 

i=l,...,N, j>i 
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Step 2. Consider now a left continuous BV function w = w(t) and let y : [0, T] ^ O be the solution 
of the corresponding Cauchy problem (1.7)-(1.8), with y G K. Observe that, for any i = 1, . . . ,N, 
and for any interval J C [0, T] such that 

y(t) 6fi\ft; for all t G J, 

the composed map <p+ o y : J R is also a left continuous BV function whose distributional 
derivative in = D (<p+ o y) is a Radon measure, which can be decomposed into an absolutely 
continuous fif c and a singular part n\, w. r. t. the Lebesgue measure dt. One can easily verify that, 
for any Borel set E C J, the absolutely continuous part of fa is given by 

H1 C {E) = f (V<p ai (y(t)), 9(y(t)) + w(t)) dt, E+ = {teE; y(t) G n at } . (2.15) 

Je+ 

Moreover, calling /x° c , fj,^, respectively the absolutely continuous and the singular part of fi w = w, 
the following bounds hold 



(V<PaMt)), «»(*)) dt + tf(E) 

E+ 



<c"'-{\n a w c (E)\ + \^ w (E)\} 

< c'" • Tot.Var.{w}, (2.16) 



for some constant c'" > that depends only on the compact set K' and on the time interval [0,T]. 
Let d, £i, i — N, N — 1 , . . . , 1 , be the constants recursively defined by 

C N = l + c"', iN = ^T> ( 2 - 1? ) 

N N 

C i =c"-£ i+1 + ]T Cj, h = -(2Ci + c"- £ o) if t<N - ( 2 - 18 ) 

Lemma 2.3. Assume that 

Tot. Var.{w} <S= (2.19) 
2Ci 

and assume that there exists some interval [t\, t 2 ] C [0, T] and some index i G {1, . . . , N} such that 

meas{t G [h, t 2 ] : a*(y(t)) = aj} < £j ■ Tot.Var.{w} V j > i (2.20)* 
together with one of the following two conditions 

(ai) 

<p ai {y{t)) < 2d ■ Tot.Var.{w} V i G [ti, t 2 ], (2.21), 

meas{t G t 2 ] : a*(y(t)) = a*} > ^ • Tot. Var.{io}. (2.22)* 



(bi) T/iere exists r G [ti, £2] sucA f/iat 

Va 4 (2/(T)) > 2 C» • To£.Far.{w}. 



(2.23), 
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Then one has 



<Pa t (v(t2)) > C\ ■ Tot.Var.{w}. 



(2.24), 



Towards a proof of the lemma, observe first that the recursive definition (2.17)-(2.18) of the 
constants Cj, li, and the bound (2.19) clearly imply 



c, > i+ c '" + c"- J2 ^ 

j=i+i 

2d ■ Tot.Var.{w} < p. 



(2.25) 



(2.26) 



Assume now that (2.20), - (2.22), hold. Then, using (2.13)-(2.16) and recalling (2.25)-(2.26), we 
obtain 



{t£[t u t 2 ] : a*(t)= ai } 



<0/(* 2 )) > + J (Vp Q , (»(*)), 9aMt))) dt " 

(V<p ai (y(t)), g aj {y(t)))\ dt ~ c "' ■ Tot.Var.{w} 

N 

J -(n ai (7r ai (y(i))), <&,,(»(*))) * - (c" ■ ]T ^ + c "') ' Tot.Var.{ W } 



S / 



j=»+l 



{t£[tl, t 2 ] : a*(t)=Q i } 
AT 

> -d-d'- e i - c "') ' Tot.Var.{w} 

j=i+i 

> C, • Tot.Var.jw}, 

proving (2.24)$. Next, assume that (2.20)$ and (2.23)$ hold, and let 

t' = sup{i G [t lt t 2 ] ■ <p at (y(t)) > 2d ■ Tot.Var.{w}}. 



(2.27) 



(2.28) 



Clearly, the bound (2.24), is satisfied if t' = t 2 since the map ip ai is left continuous. Next, consider 
the case r' < t 2 - By similar computations as in (2.27), using (2.13)-(2.16) and thanks to (2.20)$, 
(2.25)-(2.26), we get 



dt- 



{te[r',t 2 ] : a'(t)=a z } 



N 



J |(Vp Q4 (»(*)), ffa,(yW))| dt - d" • Tot.Var.{ W } 



-E 

{te[r',t 2 ] : a*(t)=a 4 } 
iV 

> (2C, - 1 - d" - c" • ^ ^) • Tot.Var.{w} 



i=»+i 



> C, • Tot.Var.-jw}, 



(2.29) 
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thus concluding the proof of Lemma 2.3. 

Step 3. Assume that the bound (2.19) on the total variation of w = w(t) holds. Set T\ = 0, 
t n+1 = T, and define recursively the points tat, tjv-i, . . . ,T2, by setting 

n = inf {i G [0, r i+ i] : <^ ai (j/(s)) > -Tot.Var.{ W } Vs £ [t, t, +1 ]}, l<i<7V. 

(2.30), 

Applying Lemma 2.3 and proceeding by backward induction on i = N, N — 1, . . . , 2, we will show 
that, for any t < r i} i = 2, . . . , N, one has 



measjs G [0, t] : a*(y(s)) — aA < £ t ■ Tot-Var-jw}, 
<P<xi(y(t)) < 2C i ■ Tot.Var.-jw}. 



(2.31), 



Indeed, if (2.31)jv is not satisfied, one of the two conditions (a) N or (b) N must be true on some 
interval [0, t], t < tn- But then, by (2.24)at, we have 

^ N (»(«)) > C N ■ Tot.Var.{w;} V a G [t, T], 

which contradicts the definition (2.30)jv- On the other hand, if we assume that (2.31)j holds for 
j = i + 1, . . . , N, but not for j = i, then one of the two conditions (a) ; or (b) ; must be true on 
some interval [0, t], t < n. Moreover, the inductive assumptions (2.31)j, j > i, imply (2.20)j and 
hence, as above, thanks to (2.24), we get 

<P ai (y(s)) > Q ■ Tot.Var.{w} Vse[t, T], 

reaching a contradiction with the definition (2.30)j. 

To conclude the proof of (i), observe that, thanks to (2.31),, i = 2, . . . , N, we have 

measjs G [n, T i+1 ] : a*{y(s)) > a,} < (X]^') ' Tot.Var.{w} V i > 1. (2.32) 

j>i 

Therefore, recalling the definitions of the map <p ai at (2.10), taking S as in (2.19), and 

N 

C>(N + l)-^2£j, (2.33) 

from (2.31)i and (2.32) we deduce that the partition n = < t 2 < • • • < r^+i = T of [0, T], 
defined at (2.30)i, satisfies the properties (2.5)-(2.7). 



2. Concerning (ii), let C, S > be the constants defined according to (i) and, given a BV function 
w = w(t) with Tot.Var.jw} < 5, and a solution y : [0, T] ^ il of the Cauchy problem (1.7)-(1.8) 
with y E K, consider the partition = t\ < T2 < ■ ■ ■ < tn + i = T, of [0, T], with the properties 
in (i). Setting 

t[ = inf {t G [n, T i+ i] : a*{y{t))=a^ i = l,...,N, 
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r(t) = 



(2.34) 



define the map 

T[ if t€]Ti,Ti] 

sup{s G [r', i] : a*(y(s)) = a,} if t g]t/, r i+ i], 

over any interval ]rj, Tj+i], i = 1,...,N. Notice that, in the particular case where a*(y(t)) > on 
for all t €]Tj, Ti+i], by the above definitions one has r(i) = r t ' = Tj+i for any t g]tj, Tj+i]. Then, let 
yO : [0, T] h-> Q be the map recursively defined by setting 



y°(*) = y(*) 



V t g]^, T], 



(2.35) 



tf°W= if Ti<Ti+i, a*(y(T(t)))>ai, V t e]n, r i+1 ], i < iV, 

< y(T(t)) if ^<n+i, a*(y(r(i))) = a;, 

2/0(0=2/0(0+) 
and let u;0 = w ^ (£) be the function defined as 

w*{t) = y<>(t) - /%(«<>(*)) ds VtG [0, T]. 
Jo 

Clearly t/O, are both BV functions as well as y, w. Moreover, yO is a solution of the perturbed 
equation (2.8). By construction, for every 1 < i < N there holds 

{ on if t[ < r i+ i, 

a*(yO(f)) = J V i efc, n +i]. (2.39) 

U*(yO(r +1 +)) if T> = n+l 

Hence the map t i— > a*(yO(^)) i s non-increasing and left-continuous. Next, recalling (2.6) and 
observing that 



(2.36) 
(2.37) 

(2.38) 



a*{y{t)) = a, 



r(t) = t, 

vHt) = y(t), 



V t g]t;, r i+ i], 



defining 



we have 



X = \J{te}n,T l+1 } : a*(y(t))> ai }, 



y(t)=yO(t) vte(o,r)\j. 

On the other hand, by the above definitions, calling M = sup ygn |g(y)|, we derive 
W{t) -t\< mcas(Z) Vtel, 



y^(t)-y(t)\< f \g(y(a))\ ds + Tot. Var.{w ; [0,t]} 

Jr(t) 



r(t) 

< M • meas(X) + Tot.Var.{w} 



Vtel, 



(2.40) 

(2.41) 
(2.42) 
(2.43) 

(2.44) 
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and 

Tbt.Var.-fj/ } - Tot.Var.{y} < Tot.Var.jy ; 1} 

< M ■ meas(J) + Tot.Var.{w}. (2.45) 

Then, using (2.44)-(2.45), we obtain 

Tot.Var.{w^} - Tot.Var.{w} < J \g(y <> (s)) \ - \g(y(s))\ ds + Tot.Var.{?/*} - Tot.Var.{y} 



< M' ■ |mcas(J) + Tot.Var.{tu}}, 



(2.46) 



for some constant M' > 0, depending only on the field g. Hence, from (2.42), (2.44), (2.46), and 
applying (2.7), it follows that y^{-) satisfies the estimates in (2.9), for some constant C > 0, which 
concludes the proof of (ii). 

We can now take 5 as in (2.19) and choose C > C according to (2.33). Both properties (i) and 
(ii) are then satisfied, completing the proof of Proposition 2.2. □ 



Proof of Theorem 1. 

For a given sequence of solutions y v : [0,T] h SI of the perturbed system (1.11) with 
Tot.Var.juv} < 8 V , b v — > as v — > oo, assume that the y v {-) converge to a function y : [0,T] 
uniformly on [0,T], and that y v (0) belongs to some compact set K C for everyy v. Thanks to 
property (ii) of Proposition 2.2, in connection with any pair w v (-), y v {-), there will be a BV function 
w$(-) and a solution y$(-) of (2.8) that satisfy 

Tbt.Var.{«£} < C ■ 5 V , \\y<> - v4^^. T] ) < C ■ 5 V , (2.47) 

for some constant C > 0, independent of v. Moreover there exists a partition = t x v < r 2 u < 

■■■< T N + 1.„ = T ° f [0) T L SUCh that 

a*(yS(t))=a i Vte]^,r j+ J, i = l,...,N. (2.48) 

Recalling (1.4) and (1.9), because of (2.48) we have 

r T e+i v r l 
y$(t) = y$(0) + J2 ' ds + / 9ai (yS(s)) ds + [w${t) - w$(0)] 

VtGlwJ, i = l,...,N. (2.49) 

By possibly taking a subsequence, we can assume that every sequence (r ; >1 converges to some 
limit point, say 

fi = lim T i v i = 1, . . . , N + 1. 



Flow Stability of Patchy Vector Fields 



11 



We now observe that 

oo oo 

m,r m [C |J f)]r it „, r i+1 J Vi. 

(1=1 V — fl 

Moreover, (2.47)2 and the uniform convergence y u (-) —> y(-) yield 

Jim ||^- y || LOO([0T]) =0. (2.50) 

From (2.47)i , (2.48)-(2.59) we now deduce 

y(t) en ai \Ufy,, 

t Vte]ri,r i+ i], Vi (2.51) 

y(t) = 2/(0) +J2 ds + /&,(»(«)) ds 

In particular, on each interval [7%, Tj+i], the function y(-) is a classical solution of y = g ai {y) and 
satisfies 

*K S ~) = .9a; V S 6]?i, T i+ l]' 

Moreover observe that, because of the inward-pointing condition (1.1), the set {i € [rj, Tj+i] : 6 
<90 ai } is nowhere dense in [r,, Tj+i]. Thus, if s is any point in ]rj, Tj+i] such that y(s) e <9f2 Qi , 
there will be some increasing sequence (s n ) n C]fi, Tj+i[ converging to s and such that y(s n ) € £l ai 
for any n. But this yields a contradiction with (1.1), because 

< hm ( VW -VM , n ai (y(s))) = (y(s-), n(y(s))) = (g ai (y('))> "«,(»(«))>■ 
Hence, recalling the definition (1.2), from (2.51) we conclude 

y(t) G Qi \|J n /5 Vte]T„f, +1 ], i = l,...,JV, 

/3>a, 

y(t) = »(0)+ / «/(»(«)) Vie [0, T], 

Jo 

proving that y : [0, T] O is a Caratheodory solution of (1.5) on [0, T\. □ 



Proof of Corollary 1.3. 

Assuming that statement is false, we shall reach a contradiction. Fix any closed subset A C O, 
any compact set K <Z A, and assume that, for some T, e > 0, there exists a sequence of solutions 
y v : [0,T] ^4 of the perturbed system (1.7), with y„(0) E K, Tot.Var.{uv} < 5 V , 5 V —> as 
v — > oo, such that the following property holds. 
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(P) Every solution x : [0,T] i— » O of the unperturbed equation (1.5) satisfies 



L°°([0,T]) 



> e 



V ZA 



(2.52) 



For each v, call y£ : [0, T] 
defined by setting 



l n the polygonal curve with vertices at the points y^(£S u ), £ > 0, 



y?(t) = y,(M„) + • (y v {{£ + l)5 v ) y v (£5 v )) 

V t E \IS V , {£ + 1)6 V ] n [0, T], < £ < [T/6 V \ , 



(2.53) 



where \T/8 V \ denotes the integer part of T/8 V . Since every ?/„(•) is a BV function that solves the 
equation (1.7), it follows that there will be some constant C > 0, independent on v, such that 



Tot.Var.^ ; J} < C ■ mcas(J) + Tot.Var.{w„ ; J} 



(2.54) 



for any interval J C [0, T]. Then, using (2.54), we derive for any fixed < I < £' < [T/5 V \ the 
bound 

V^S„) - y$(t5 v )\ = \y v (t'5 v ) - y v {td v ) 

< Tot. VaT.{y u ; [£5 V , £'8 V ]} 

< (1 + C) ■ {I' - l)8 v . (2.55) 

Therefore y$(-) is a uniformly bounded sequence of Lipschitz maps, having Lipschitz constant 
Lip(j/5 > ) < (1 + C). Hence, applying Ascoli-Arzela Theorem, we can find a subsequence, that we still 
denote y$(-), which converges to some function y : [0, T] M™, uniformly on [0, T]. On the other 
hand, by construction and thanks to (2.54), for any fixed < t < T, with £S U <t< (£+ 1)5 V , there 
holds 

y„(t)-y$(t)\ < \y v {t) - y v {£5 v )\ + \y v {£8 v ) - y${t) 

< y v {t) - y v {iS v )\ + \y v ({£+l)5 v ) - y u {£5 v ) 

< 2-Tot.Var.{j/, ; [£S V , (t+l)5 v ]} 

< 2(1 + C)-5 V . (2.56) 
Thus, since 8 V — > as v — > oo, the uniform convergence of y„(-) to y(-) implies 

lim \\y v - y 



lLoo([ ,T]) 



0. 



(2.57) 



By assumption, Range(y J/ ) C AcO for every v, and hence from (2.57) we deduce that also the limit 
function y(-) takes values inside fi. We can thus apply Theorem 1 to the sequence ?/„(•) and conclude 
that the function y : [0, T] il is a Caratheodory solution of the unperturbed equation (1.5) with 



\\y 2mIl°°([o,t]) < 6 

for all v sufficiently large. We thus obtain a contradiction with (2.52), concluding the proof. 



□ 
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3 - Robustness of Patchy Feedbacks. 

In this section we apply the previous results on patchy vector fields with impulsive perturbations, 
and construct (discontinuous) stabilizing feedback controls that enjoy robustness properties in the 
presence of measurement errors and external disturbances. Consider the nonlinear control system 
on R" 

x = f(x,u) u(t) G /C, (3-1) 

assuming that the control set /C C W n is compact and that the map / : R™ x M. m >—> R™ is smooth. 
We seek a feedback control u = U{x) e JC that stabilizes the trajectories of the closed-loop system 

x = f(x, U(x)) (3.2) 

at the origin. It is well known that, even if every initial state x £ R™ can be steered to the origin by 
an open-loop control u = u x (t) , a topological obstruction can prevent the existence of a continuous 
feedback control u = U(x) which (locally) stabilizes the system (3.1). This fact was first pointed out 
by Sussmann [Su] for a two-dimensional system (n = 2, JC — R 2 ), and by Sontag and Sussmann [SS] 
for one-dimensional systems (n = 1, K = R). For general nonlinear systems, it was further analyzed 
by Brockett [Bro] and Coron [Corl]. It is thus natural to look for a stabilizing control within a class 
of discontinuous functions. However, this leads to a theoretical difficulty, because, when the function 
U is discontinuous, the differential equation (3.2) may not admit any Caratheodory solution. To 
cope with this problem, two different approaches have been pursued. 

1. An algorithm is defined, which constructs approximate trajectories in connection with an arbi- 
trary (discontinuous) feedback control function. For example, one can sample the feedback control 
at a discrete set of times. In this case, one is not concerned with the existence of exact solutions, 
but only in the asymptotic stabilization properties of all approximate solutions. 

2. Alternatively, by the asymptotic controllability to the origin of system (3.1) by means of open- 
loop controls, one proves the existence of a stabilizing feedback u = U (x) having only a particular 
type of discontinuities. This feedback thus generates a patchy vector field, and the corresponding 
system (3.2) always admits Caratheodory solutions. 

The first approach was initiated in [C-L-S-S], and further developed in [Ril], [Ri2]. The second 
was introduced in [A-B], defining the following class of piecewise constant feedback controls: 

Definition 3.1. Let (O, g, (Q a , g a ) a€A ) be a patchy vector field. Assume that there exist control 
values k a e 1C such that, for each a e A, there holds 

9a{x) = f(x, k a ) V X G D a = ft Q \ (J Sip. (3.3) 

I3>a 

Then, the piecewise constant map 

U(x) = k a if x e D a (3.4) 
is called a patchy feedback control on 0, and referred to as (ft, U, (fi Q , k a ) aeA ) . 
Remark 3.2. From Definitions 1.2 and 3.1, it is clear that the field 



g(x) = f(x, U(x)) 
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defined in connection with a given patchy feedback (fi, U, (fl a , fc Q ) aeA ) is precisely the patchy 
vector field (O, 3, (fi a , <7a) a€A ) associated with a family of fields {g a : a <E A) satisfying (1.1) 
Clearly, the patches (fi Q , g a ) are not uniquely determined by the patchy feedback U. Indeed, 
whenever a < (3, by (3.3) the values of g a on the set Q a \ are irrelevant. Moreover, recalling the 
notation (1.3) we have 

U(x)=k a .( x) Vieil. (3.5) 

Here, we address the issue of robustness of a stabilizing feedback law u — U(x) w. r. t. small 
internal and external perturbations 

x = f(x, U(x + ((t)))+d(t), (3.6) 

where £ = £(t) represents a state measurement error, and d = d(t) represents an external disturbance 
of the system dynamics (3.2). Since we are dealing with a discontinuous O.D.E., one cannot expect 
the full robustness of the feedback U(x) with respect to measurement errors because of possible 
chattering behaviour that may arise at discontinuity points (see [Hel], [So4]). Therefore, we shall 
consider state measurement errors which are small in BV norm, avoiding such phenomena. 

Before stating our main result in this direction, we recall here some basic definitions and Propo- 
sition 4.2 in [A-B] . This provides the semi-global practical stabilization (steering all states from a 
given compact set of initial data into a prescribed neighborhood of zero) of an asymptotycally 
controllable system, by means of a patchy feedback control which is robust with respect to ex- 
ternal disturbances. We consider as (open- loop) admissible controls all the measurable functions 
u : [0, 00) -> R m such that u(t) E JC for a.c. t > 0. 

Definition 3.3. The system (3.1) is globally asymptotically controllable to the origin if the following 
holds. 

1. Attractiveness: for each x € R™ there exists some admissible (open-loop) control u = u x (t) 
such that the corresponding trajectory of 

x(t) = f(x(t),u*(t)), x(0)=x, (3.7) 

either reaches the origin in finite time, or tends to the origin as t —> 00. 

2. Lyapunov stability: for each e > there exists 6 > such that the following holds. For 
every x G K™ with \x\ < 5, there is an admissible control u x as in 1. steering the system from 
x to the origin, so that the corresponding trajectory of (3.7) satisfies |x(£)| < e for all t > 0. 

Proposition 3.4. [A-B, Proposition 4.1] Let system (3.1) be globally asymptotically controllable to 
the origin. Then, for every < r < s, one can find T > 0, X > 0, and a patchy feedback control 
U : D 1 — > /C, defined on some domain 

D D {x e R" ; r < \x\ < s} (3.8) 

so that the following holds. For any measurable map d : [0, T] \— ► R" such that 

IHIl°°([o,t]) - * ! 
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and for any initial state xq with r < \xq\ < s, the perturbed system 

x = f(x,U(x)) +d(t) (3.9) 

admits a (forward) Caratheodory trajectory starting from x - Moreover, for any such trajectory 
t i ► 7(i), t > 0, one has 

7(i) e D V i > 0, (3.10) 



and £/iere existe t 7 < T such that 



| 7 (t 7 )|<r. (3.11) 



Relying on Corollary 1.3 of Theorem 1 and on Proposition 3.4., we obtain here the following 
result concerning robustness of a stabilizing feedback w. r. t. both internal and external perturba- 
tions. 

Theorem 2. Let system (3.1) be globally asymptotically controllable to the origin. Then, for every 
< r < s, one can find T" > 0, x' > 0, and a patchy feedback control U' : D 1 1— ► K defined on some 
domain D' satisfying (3.8), so that the following holds. Given any pair of maps ( S BV([0, T']), 
d e L°°([0, T']), such that 

Tot.Var.{(} < X ' , IMIl~([o,t<]) ^ x ' ' ( 3 ' 12 ) 

and any initial state xq with r < \xq\ < s, for every solution t ^ x(t), t > 0, of the perturbed 
system (3.6) starting from Xo, one has 

x{t) efl' V t e [0, T'] , (3.13) 

and there exists t x <T', such that 

|x(*x)|<r. (3.14) 



Proof. 

1. Fix < r < s. Then, according with Proposition 3.4, we can find T' > 0, and a patchy feedback 
control U' : D' 1— > K,, defined on some domain 

D' D {x e K™ ; r/3 < \x\ < s} , (3.15) 

so that the following holds. For every Caratheodory solution t >— > x(t), t > of the unperturbed 
system (3.2) (with U = U' ) starting from a point x in the compact set 

K = {x e K" ; r < \x\ < s}, (3.16) 

one has 

x(t) eD p = {xeD' : d(x, dD') > p) V t > 0, (3.17) 
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for some constant p > 0. Moreover, there exists t x < T' such that 

\x(t x )\ < (3.18) 

According with Definition 3.1, the field 

g(x) = f(x, U'(x)) (3.19) 

is a patchy vector field associated to the family of fields {g a : a £ A] defined as in (3.3). The 
smoothness of / guarantees that, for BV perturbations w — w(t) having some uniform bound 

Tot.Var.jw} < £ on the total variation, every (left-continuous) solution y : [0, T'] i— > R 2 of the 
impulsive equation (1.7), starting at a point x 6 K, takes values in the closed set 

A = B(D p , p/2). (3.20) 

Therefore, thanks to Corollary 1.3 of Theorem 1, there exists some constant 

< x' = X~'(A K, r, r/3) < x (3-21) 

such that the following holds. If y : [0,T'] R 2 is a (left-continuous) solution of the impulsive 
equation (1.7), with y(0) e K and Tot.Var.(w) < then one has 

y(t) e A V t e [0, T'] , (3.22) 

and there exists t y < T' such that 

\y(t y )\ < y- (3.23) 

2. In connection with the patchy feedback V introduced above, define the map 

h(y, z) = f(y - z, U'{y)) - f(y, U'(y)) (3.24) 
and observe that, by the smoothness of /, there will be some constant c > 0, such that 

\h(y,z)\ <c-\z\ VyeA, \z\ < x' ■ (3.25) 
Consider now a pair of maps ( G BV([0, T']), d e L°°([0, T']), satisfying (3.12) with 



X<™in\ ^ X , v I, P -\, (3.26) 



2(1 + T'c'Y 3' 2 

and let x = x(t) be any Caratheodory solution of the perturbed system (3.6), with an initial 
condition x(0) = x £ K. Then, as observed in the introduction, the map 

t h+ y(t) = x(t) + C(t) (3.27) 
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satisfies the impulsive equation (1.7) where 

w(t) = C(*) + / (%(*), C(s)) + d(s)) ds . (3.28) 
Jo 

But then, since (3.12), (3.25), (3.26), together, imply 

Tot.Var.{ W ; [0, T']} < Tot.Var.jC ; [0, T'}} + T'c ■ \\C\\^ {[0 TI]) + T' ■ \\d\\^ {[0 TI]) 
< (1 + T'c) • Tot.Var.jC ; [0, T']} + \\d\\^ {[a TI]) 

<x', 

from (3.22)-(3.23) and (3.12), (3.20), (3.26), (3.27) it follows 

x{t) G B(i4, p/2) CD', V t g [0, T'] , 
| a;( t y ) | < r for some t v <T', 
which completes the proof of the theorem, taking \' as m (3.26). □ 



(3.29) 



Remark 3.5. For discontinuous stabilizing feedbacks constructed in terms of sampling solutions, 
an alternative concept of robustness was introduced in [C-L-R-S], [So4]. In this case, one considers a 
partition of the time interval and applies a constant control between two consecutive sampling times. 
To preserve stability, the measurement error should be sufficiently small compared to the maximum 
step size. Moreover, each step size should be big enough to prevent possible chattering phenomena. 
The next result shows that the feedback provided by [A-B, Proposition 4.2] enjoys also this type of 
robustness. Before stating this result we describe now the concept of sampling trajectory associated 
to the perturbed system (3.6) that was introduced in [C-L-R-S]. 

Let an initial condition xq, and a partition tt = {0 = To < t± < ■ ■ ■ < r m+ i = T} of the 
interval [0, T], be given A sampling trajectory x w of the perturbed system (3.6), corresponding to 
a set of measurement errors {ei}™ and an external disturbance d G L°°([0, T]), is defined in a 
step-by-step fashion as follows. Between To and n, let x % (-) be a Caratheodory solution of 

x = f{x, U(x + e )) + d(t) t G [to, n] , (3.30) 

with initial condition ^(O) = x a . Then, x n (-) is recursively obtained by solving the system 

x = f(x,U(x*(Ti) + ei))+d(t) te[n, Ti+i], i>0. (3.31) 

The sequence {x^^i) + ei}™ corresponds to the non-exact measurements used to select control 
values. 



Theorem 3 Let system (3.1) be globally asymptotically controllable to the origin. Then, for every 
< r < s, one can find T" > 0, x" > 0, <5 > 0, k > 0, and a patchy feedback control U" : D" i— > K, 
defined on some domain D" satisfying (3.8) so that the following holds. Given an initial state xq 
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with r < \xq\ < s, a partition ir = {tq = 0, t\, ... , T m +i = T"} of the interval [0, T"} having the 
property 

2<T i+ i—Ti<5 Mi, for some 5 S ]0, 5], (3.32) 

a set of measurement errors {ei}™ and an external disturbance d G L°°([0, T"]) that satisfy 

max \ei\ < k ■ 5 , (3.33) 

i 

IMI L oo<x", (3-34) 
the resulting sampling solution x 7r (-) starting from x has the property 

x v (t)£D" V t £ [0, T" ] . (3.35) 
Moreover, there exists t x ^ < T" , such that 

\x n (t x M<r. (3.36) 



Proof. 

1. Fix < r < s. Then, according with Proposition 3.4, we can find T' > 0, \' > 0, and a patchy 
feedback control U" : D" i — > K, defined on a domain 

D" R" ; r/3 < \x\ < 2s} 

so that the following holds. For every external disturbance d G L°° satisfying (3.34) with \" < x' , 
and for any Caratheodory solution t ^ x(t), t > of the perturbed system (3.9) (with U = U"), 
starting from a point x with r < \x n \ < s, one has 

x{t) e D P1 = {x e D" : d(x, dD") > Pl } V t > 0, (3.37) 

for some constant p\ > 0. Moreover, there exists t x < T' such that 

\x(t x )\ < T -. (3.38) 

Let 

) : a = l,...,N} g a (x) = f(x,k a ), k a G K, (3.39) 

be the collection of patches associated with the patchy vector field 

g(x) = f(x,U"(x)). (3.40) 

We may assume that every vector field g a is defined on a neighborhood B(Q a , p 2 ), < p2 < pi of 
the domain fi Q so that, setting 

n p a = {xen a ; d(x, on a ) > P ) , (3.41) 
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one has 

and that every g a is uniformly non-zero on the domain D a defined in (3.3). Moreover, thanks to 
the inward-pointing condition (1.1), we may choose the constants < p 2 < r/3, and x" < x' so 
that there holds 

|<7a(a;)| >2 X " VxeB(D a ,p 2 ), (3.42) 

and 

(g a (x) + v, n(aO) < Vxe B(dQ a , p 2 ), \v\ < X " • (3.43) 

For every d G L°°, we denote by t h-> x a (t; to, x , d) the solution of the Cauchy problem 

x = g a (x) + d(t), x(t ) = x , (3.44) 

and let [to, i max ] be the domain of definition of the maximal (forward) solution of (3.45) that is 
contained in B(D a , p 2 ). 

Observe that, since every Caratheodory solution of the perturbed system (3.9) (with U = U"), 

o 

starting from a point xq e -8(0, s)\B(0, r), reaches the interior of the ball B(0, r/3) in finite time, 
and because of (3.42), for any a = 1, . . . , N one can find T a > with the following property. 

(P)i For every x E B(D a , p/2), < p < p 2 , and for any d G L°° satisfying (3.34), there exists 
some time t p = t p (xo, d) < T a such that, either one has 

2r 

\x a (t + t p ;t , x , d) \ < —, (3.45) 

or else there holds 

x a (t;t , x , d) eB(D a , p 2 )\B(D a , p) V t G [t + t p , i max ]. (3.46) 

On the other hand, relying on the inward-pointing condition (3.43), we deduce two further properties 
of the solutions of (3.44). 

(P)2 The sets < p < p 2 , defined in (3.41) are positive invariant regions for trajectories 
of (3.44), i.e., for every xo G and for any d G L°° satisfying (3.34), one has 

x a (t- t , x , d) G W a Vt>i - (3.47) 

(P)3 There exists some constant c > so that, for every xq G B(il a , p), < p < p 2 , such that 
d(xo, dSl a ) < p, and for any d G L°° satisfying (3.34), one has 

x a (t; t , x , d) G n 2 a p Vt>t + c-p. (3.48) 



2. Consider an initial state xo G B(0, s)\B(0, r), and a partition ir = {nji^o OI [0, oo[ having 
the property (3.32), with 

0<5<5 = mm\c-p 2 , g}, (3.49) 
M = sup {\g a (x)\ : x £ B(n a , p 2 ), a = l,...,N}. 
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Let Xt, : [0, oo[ i— > M™ be a sampling solution starting from xq, and corresponding to a set of 

measurement errors {ei}™ and to an external disturbance d(-) G L°° that satisfy (3.33)-(3.34) 
with 

fc=^. (3.50) 

We will first show the following 
Lemma 3.6. The map 

ii — >a*{Ti) = a*(x v {Ti) + ei) z>0, (3.51) 

is non- decreasing. 

Indeed, assume that a*(rj) = a, which, by definitions (1.3), (3.3), (3.5) implies 

x n (n) + e l G Da , (3.52) 

x n {n+i) = x a (r i+ i] n, x n (Ti), d t[ TilTi+1 ] ) , (3.53) 
Then, because of (3.33), (3.49)-(3.50), one has 

Xi = x n (Ti) e B(D S , kg) c B(fi a , pa). (3.54) 
We shall consider separately the case in which 

x t G Df c nf kS< P2 , (3.55) 

and the case where 

x t G B{D dl , kS), d(xi , dQ & ) < kS < p 2 . (3.56) 

In the first case, using (3.53) and applying (P)2 we deduce that x n (Ti+i) G 0~ 5 which, in turn, 
because of (3.33), (3.49)-(3.50), implies 

XTr(T i+ i) + ej+i G fia- (3.57) 

From (3.57), by definition (1.3) we derive 

a*(T»+i) > S , (3.58) 

proving the lemma whenever (3.55) holds. On the other hand, when (3.56) is verified, since by 
(3.32), (3.50) one has 

applying (P)3 we deduce 2^(7^+1) G fi a fe5 . This again implies (3.57)-(3.58), completing the proof 
of Lemma 5.6. 

Next, relying on (P)i, and setting 

i a = min{i>0 ; a*(r 4 )=a, x n (n) B(0, 2r/3) } , 

a G Range(a*), 

«'^=max{i>0; a*(n)=a, x^n) <£ 5(0, 2r/3) } , 

(3.59) 
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we deduce 

t.„-t., <T a V a G Range(a*) . (3.60) 

Indeed, if (3.60) does not hold, by definitions (3.3), (3.5) one has 

xv = x„(t, ) G B(D a , kS) c B(fi Q , pa/2) , (3.61) 
x w (t) = x a (t; T t , a ,x Va , d \ [T{i ) , V t G [r^, r 4 „ + J . (3.62) 

But then, applying (P)i, one could find some i < i" a such that 

x v (t) G B(D a , p 2 ) \ B{D a , 2kS) V t G [r ? , t ( „ +1 ] . 

By definitions (1.3), (3.51) and because of (3.33), this implies 

a*(n)>a Vi<i<il, 

providing a contradiction with (3.59). 

To conclude the proof of Theorem 3, we observe that the monotonicity of the map (3.51), 
together with the estimate (3.60), implies that there exists some time t x „ < T" = Yl2=i T « such 
that (3.36) is verified. Moreover, (3.35) clearly follows from (3.37) and (3.49). □ 

Remark 3.7. Consider a partition tt = {r = 0, n, ... , r TO+ i = T} of the interval [0, T] having 
the property (3.32). If we associate to a set of measurement errors {e^}™! satisfying (3.33) the 
piecewise constant function ( : [0, T] i— > W 1 defined as 

C(i) = e { Vt G ]n, T i+ i] , 

then 

Tot.Var.{C} < Ak ■ T . 

Thus, taking the constant k sufficiently small we may reinterpret the discrete internal disturbance 
allowed for a sampling solution in Theorem 3 as a particular case of the measurement errors with 
small total variation considered in Theorem 2. 
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